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Abstract 

We present a novel set of Feynman rules and generalised unitarity cut-conditions for computing one-loop amplitudes 
via c/-dimensional integrand reduction algorithm. Our algorithm is suited for analytic as well as numerical result, 
because all ingredients turn out to have a four-dimensional representation. We will apply this formalism to NLO QCD 
corrections. 
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1. Introduction 

In the era of the LHC experiments of increasing ac¬ 
curacy become possible, where one of the highlight of 
Run 1 of the LHC was the discovery by CMS and AT¬ 
LAS of a Higgs boson dll. Hence it is necessary to 
achieve more accurate results for measurable quantities 
at the theoretical level. 

According to perturbation theory, higher order correc¬ 
tions to amplitudes have to be considered. To evaluate 
such corrections in quantum field theory, it is necessary 
to compute multi loop Feynman diagrams, where, in¬ 
stead of the explicit set of loop Feynman diagrams, the 
basic reference point is the linear expansion of the am¬ 
plitude function in a basis of master integrals (Mi’s), 
multiplied by coefficients that are rational functions of 
the kinematic variables, already known as Passarino 
Veltman reduction theorem d at the level of one-loop. 
It is in fact possible to recover the finiteness of scattering 
amplitudes at integrand level by constructing the inte¬ 
grands by a multi-particle pole expansion arising from 
the analiticity properties and unitarity of the S-matrix. 
Indeed, scattering amplitudes, continued for complex 
momenta, reveal their singularities structures as poles 
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and branch cuts. The unitarity based method (UBM) 
allows to determine the coefficients of the Mi’s by ex¬ 
panding the integrand of the tree level cut amplitudes 
into an expression that resembles the cut of the basis 
integrals. 

In this talk, I review the four dimensional formula¬ 
tion (FDF) proposed in (01 which is equivalent to the 
four-dimensional helicity (FDH) scheme Isl-Ql, and al¬ 
lows for a purely four-dimensional regularisation of the 
amplitudes. Within FDF, the states in the loop are de¬ 
scribed as four dimensional massive particles. The four¬ 
dimensional degrees of freedom of the gauge bosons are 
carried by massive vector bosons of mass /r and their 
(d - 4)-dimensional ones by real scalar particles obey¬ 
ing a simple set of four-dimensional Feynman rules. A 
c/-dimensional fermion of mass m is instead traded for 
a tardyonic Dirac field with mass m + ipy^ |0]. The 
d dimensional algebraic manipulations are replaced by 
four-dimensional ones complemented by a set of mul¬ 
tiplicative selection rules. The latter are treated as an 
algebra describing internal symmetries. 

This contribution is organised as follows; sectionl^is 
devoted to the description of the regularisation method, 
while Section |3] describes how generalised unitarity 
method can be applied in presence of a FDF of one- 
loop amplitudes. Section |4] shows the decomposition in 
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terms of Mis of certain classes of 2 —> (n — 2) one-loop 
amplitudes. Section |5] collects the applications of gen¬ 
eralised unitarity methods within the FDF. In particular 
it presents results for representative helicity amplitudes 
of gg gg, gg ggg, gg gggg and gg gH. 

2. Four dimensional formulation of the d- 
dimensional regularisation scheme 

We discuss briefly the regularisation scheme pro¬ 
posed in 101 , only pointing out the main ingredients. 

Let’s denote as barred a quantity referred to unob¬ 
served particles and living therefore in a ^-dimensional 
space. Then, the metric tensor can be split as 

= ( 1 ) 

in terms of a four-dimensional tensor g and a - 2 e- 
dimensional one, g, such that 

r^gpy^Q, 4=4, ( 2 ) 

The tensors g and g project a c/^-dimensional vector q 
into the four-dimensional and the - 2 e-dimensional sub¬ 
spaces respectively. 


3. Generalised Unitarity 

In this section we discuss the consequence of using 
internal lines in (4 - 2e)-dimensions within FDF where 
spinors and polarisation vectors are written explicitly. 
These ingredients allow the construction of the tree- 
level amplitudes that are needed to recover any one-loop 
amplitude. 

Due to FDF scheme is suitable for the four-dimensional 
formulation of c/-dimensional generalised unitarity, all 
kinematics in the construction of the amplitude admit an 
explicit representation in terms of generalised spinors 
and polarisation expressions. 

In the following discussion we will decompose a d- 
dimensional momentum ( as follows 

+ (7) 

while its four-dimensional component { will be ex¬ 
pressed as 

2 2 

t - + q(, qt= ^ qt , ( 8 ) 

2.1-qt 

in terms of the two massless momenta C’ and q(. 


q^ = ^vq\ q^ = ^vq''- (3) 

and the properties for the matrices = g^vf'' can be 
obtained from Eq. (I2]i 

[4r']-0, {4/) = o, (4a) 

{r, 4 - 24 . (4b) 

In principle, we could infer the behaviour of 7 “ from 
dUi and say y ~ y^, however, this choice does not fulfil 
the Clifford algebra when ds 4. It means we cannot 
have any four-dimensional representation of the - 2 e- 
subspace, therefore, we introduce an algebra with an in¬ 
ternal symmetry called -2e selection rules, (-2e)-SRs, 
which consists in performing the substitutions 

4 ^ r , f (5) 


The -2e-dimensional vectorial indices are thus traded 
for (-2e)-SRs such that 


qABqBC _ 

Y'A qAB _ 


qAA ^ qAB ^ qBA^ 

FArA^O, e^F^ = l, 

= 1 . ( 6 ) 


The exclusion of the terms containing odd powers of 
ji completely defines the FDF, and allows one to build 
integrands which, upon integration, yield to the same 
result as in the FDH scheme. 


3.1. Spinors 

The spinors of a (f-dimensional fermion have to fulfil 
a completeness relation which reconstructs the numera¬ 
tor of the cut propagator, 

2 (*- 2)/2 

^ “-l,(rf) (^) (f) = / -t- ™ , 

A=l 

2 (*- 2)/2 

2 V,\, (rf) (f) V,1, (rf) (f) = / - m . (9) 

/l=l 

The substitutions (|5]l allow one to express Eq. (|9l) as 
follows: 


^ Mi {£) Ua + ipy^ + m , 

^=± 

^ Vi (f) Vi (f) = / H-///y^ - m. ( 10 ) 

^=± 


with the generalised massive spinors H 
«+(2 = 4 + ’ u -{() = \ C ] + 


(m + ip) 
(m + ip) 


\qe), 

\qt). 


(Ha) 


M+ (f) = [^\ + 


(m + ip) 
{qe^') 


{qt\ , 


u_ {()^(C\ + 


(m — ip) 


[qe\ , 
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V- {£) = [^^1 


{m + ifi) 


iqA , 


v + (i) = (^^1 


(m — iy) 


[qi\ . 


(lib) 


fulfil the completeness relation (fTOl) . 


3.2. Polarisation vectors 

In the axial gauge, the helicity sum of the transverse 
polarisation vector is 


d-2 




-r + 


P if + rtf 
T~ri 


( 12 ) 


where f) is an arbitrary ^f-dimensional massless momen¬ 
tum such that £ ■ i) 

In particular the choice 


(13) 

with (, £ defined in Eq. (I7]i, allows us to disentangle the 
four-dimensional contribution form the c/-dimensional 
one: 




(14) 


The first term is related to the cut propagator of a mas¬ 
sive gluon and can be expressed as follows 

mpv 

^ ,t =±,0 


in terms of the polarisation vectors of a vector boson of 
mass ji 








V 2 p 


(16) 


These polarisation vectors are orthonormal and display 
all of the usual properties expected for massive vector 
bosons 


The factor can be easily accounted by defining the 
cut propagator as 


(19) 

b,B ^ ' 


From generalised spinors and polarisation vectors the in¬ 
dependence of the tree-level amplitude arises. 

The PDF approach to reconstruct the rational part of 
one-loop scattering amplitudes is different from the su¬ 
persymmetric decomposition in and from the six¬ 
dimensional formalism d. Indeed, to compute any 
one-loop amplitude via supersymmetric decomposi¬ 
tion one splits the amplitude in two terms: i) cut 
constructible part which is obtained by using four¬ 
dimensional unitarity, ii) and the rational one that is 
reached by introducing in the amplitude a complex 
scalar in ^-dimensions and deal with a massive four¬ 
dimensional ones. 

On the other hand, the six-dimensional helicity method 
treats ^f-dimensional on-shell momenta into a six¬ 
dimensional massless basis and, on the cuts, uses six 
dimensional helicity spinors to compute the relevant 
tree-level amplitudes. However, because of the argu¬ 
ment given in ifl^ . the contribution that comes from this 
treatment gives a result that has to be corrected by hand 
with the help of topologies involving complex scalars 
along the lines. 

Unlike the approaches presented above, FDF does not 
make any distinction between cut-constructible or ratio¬ 
nal part, as well, the result obtained with FDF scheme 
is automated corrected by the (-2e)-SRs, it splits the 
^f-dimensional objects into their four-dimensional and 
{d - 4)-dimensional parts and finds a four-dimensional 
representation for both of them. Moreover, the ap¬ 
proaches already described are simpler than the ones 
that introduce explicit higher-dimensional extension of 
either the Dirac 1 12, lH or the spinor ifTll[T3] algebra. 


4. One-loop amplitudes 


0 , s^{€) ■ s^{€) ^, 

eoW = -l, £±(0-eo(0= 0, 

0. (17) 

The second term of the r.h.s. of Eq. (fT4l) is related to the 
numerator of cut propagator of the scalar Sg and can be 
expressed in terms of the (-2e)-SRs as: 

^ = G'^^ - . (18) 

P 


In order to apply generalised-unitarity methods 
within FDF, we consider as examples the one-loop 2 —» 
2,3,4 scattering amplitudes, where external particles 
are gluons. 

In general, due to the reduction theorem any mass¬ 
less four-point one-loop amplitude can be decomposed 
in terms Mis, as follows 

1 " * r 

‘°°^ = , ■ ' Y \ci\j\k\l-,0h\j\k\l +Cij\k\l-Qli}\k\l 
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■Cij\kl-,oIij\kl 
1 


(4;r)2 


+ 'R, ( 20 a) 

Ci\j\k\l;4li\j\k\l[f^'^] + {Cij\k\l;2 Iij\k\l[f4^] 


■ Cijlkl-llijlklir ] 


( 20 b) 


The quadruple cuts read as follows 


(-[ 0 ] _ 
^l|2|3|4 




In Eq. (l20l i. we see the decomposition between 
cut-constructible and rational part, where the latter has 
been collected in 'R. However, we emphasise one-loop 
processes are not computed by distinguishing those two 
pieces, instead within the PDF the two contributions 
are computed simultaneously from the same cuts. 


rW 

'^l|2|3|4 


^[ 2 ] 

^1|2|3|4 



(23a) 

(23b) 


Z : 



+ C.p. , 


The coefficients c’s entering in the decomposi¬ 
tions (l 20 l) can be obtained by usi^ the generalised 
unitarity techniques for quadruple |3 3 . triple ifl^ 
18], and double lfl9l - l3 cuts. Since internal particles 


are not 


are massless the single-cut techniques 11221 - 
needed. In general, the cut C,, . . 4 , defined by the condi¬ 
tions D,, = ■ ■ ■ - Dii^ - 0 , allows for the determination 
of the coefficients 


^[3] 


y T3 i i -Hc.p., 


rW 


= 7-4 



1+ 4+ 


(23c) 

(23d) 


(23e) 


where the abbreviation “c.p.” means “cyclic permuta- 
5. Examples tions of the external particles”. In Eqs. (I23t , the (-2e)- 

SR have been stripped off and collected in the prefactors 
5.7. The all plus four-gluon amplitude 7 ^^ 


First, let us consider one-loop four-point amplitudes 
with four outgoing massless particles 

0 ^ 1(Pi)2(p2)3(p3)4(p4), ( 21 ) 

where p, is the momentum of the particle i. 

Within the PDF, we consider the colour-ordered 
Feynman rules that contain interactions between gluons 
and scalars, however, due to the (-2e)-SRs, the relevant 
interactions are; /) three gluons and if) one gluon with 
two scalars, see the discussion below. 

Fet us compute the four-gluon colour-ordered helicity 
amplitude A 4 (l^, 2 ^, 3 g, 4 ^), which at tree-level van¬ 
ishes, while the one-loop contribution is finite and is 
obtained from the quadruple-cut Cippp. 

Since contribution to this amplitude comes only from 
the boxes and in FDF we have five boxes, we decom¬ 
pose this sum of boxes as: 

4 4 

Ci| 2|3|4 = ^ ^1|2|3|4 ’ ‘3l|2|3|4;n = ^ ‘'1|2|3|4 ;h ’ 

(=0 (=0 

where (c^*') is the contribution to the cut (coefficient) 
involving i internal scalars. 


T, 


= 0, 

T2 


= 0, 

7-3 


= 0, 

7-4 

= tr(GGGGGGGG) 

= - 1 


The prefactors T\,.. .,T^ force the cuts (I23bl) - (I23dl i 
to vanish identically. The only cuts contributing, 
Eqs. (I23ab and (I23eb . lead to the following coefficients 


c = 0 

‘'1|2|3|4;0 ’ 

= 0 

‘'1|2|3|4;0 ’ 


„[ 0 ] 

'1|2|3|4;4 


„[4] 

'1|2|3|4;4 


. [12] [34] 
'<12X34)’ 
. [12] [34] 
'<12X34) ■ 


(25) 


Therefore the only non-vanishing coefficient, ci( 2 | 3 | 4 ; 4 , 
is 


JOI , J41 .,,[12] [34] 

Cl|2|3|4; 4 - Cj|2|3|4. 4 + Ci|2|3|4. 4 “ ^ J2) <34) ' 

The colour-ordered one-loop amplitude can be obtained 
from Eq. (EOll, which in this simple case reduces to 

"^4(lg>2^,3j,4g) = Ci|2|314;4 h\l\-i\4\fl^] 
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i [12] [34] 
48 7:2 <12)<34> ’ 


(27) 


and is in agreement with the literature 1125 1. 

It is worth to mention the (-2e)-SRs can be performed 


before the tree-level amplitudes are computed as well 
automated cut-by-cut. It means topologies which 
vanish because of the (-2e)-SRs can be removed at 
the beginning of the computation without affecting our 
result. 


In ||3l the method is also checked by computing the 
QCD one loop correction to the scattering amplitude of 
two gluons production by quark anti-quark annihilation 
and the processes at two-loop of three gluons fusion in 
a Higgs. 

5.2. Five- and six-gluon amplitudes 

In this section we show how PDF works at higher 
point amplitudes, as an explicit example we present 
the analytic contribution of the five- and six- point all 
plus amplitude, we choose this helicity configuration 
because of the absence of triangles and bubble contri¬ 
butions. 

5.2.1. Five gluon amplitudes 


triangles and bubbles arise, then, to obtain these con¬ 
tributions we consider the topologies showed in Fig. [1] 
with the following trees needed as input: 


GG^g, 

SS^g, 


GG -> gg. 

s s gg, 


GG -> ggg. 

s s ^ ggg. 

(30) 


Where G and S are the generalised gluon and colour 
scalar respectively and, g represents the external gluon. 
As was discussed in section EU tree levels containing 
both generalised gluon and colour scalar as internal legs 
do not contribute to the coefficient due to the -2e-SRs. 
It means that to recover any five-gluon amplitude for a 
particular helicity configuration we compute the coef¬ 
ficients that appear in eq. (l20l) . obtaining an agreement 
with NJet 12711 and reproducing previous results 12811 . 



Figure 1: Triangle and bubble topologies for the five-point. 



Cl2|3|4|5;0 

Cl2|3|4|5;4 


0 , 

2/[21] [43] [53] [54] 
<12)tr5(4,l,5,3) ’ 


(28) 


with trj (1,2,3,4) = <11234| 1] - [11234| 1). 

From eq. (i20l) . the finite colour-ordered one-loop ampli¬ 
tude reduces to 


^1 '°°P(1 + ^2^,3^,4^,5^) - Ci2|3|4|5;4 7i2|3|4|5 

-I- cyclic perms. (29) 


5.2.2. Six gluon amplitudes 


Ci23|4|5|6 - 



Cl23|4|5|6;4 


2i [56] 


<12)<23)tr5 (5,4,6,l)tr5(5,4,6,3 ) 
x(s45 (6|1+2|3][51][64]2 


^46<5|1+2|3] [54]2[61] 




(32a) 


In agreement agrees with lE^ . 

It is worth to mention that within PDF we 
have also computed other helicity configurations, 
A5(r,2+,3+,4+,5+), A5(r,2-,3+,4+,5+) and 
As (1“, 2^, 3“,4^, 5^), where contributions from 
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2!<5|1 +2|6] <6|1 +2|5][12] [43] [65]^ _ 

“ <12)<23)tr5(5,2,6, l)tr5(5,4,6,3) ’ 

(32b) 



Cl2|3|45|6;0 - 0 , 

2/[12][54][63]2 

“ <12> <45> trs (2,3,6,1) trs (5,3,6,4) 
x«3|l+2|3](6|l+2|6]-536*12); 

(32c) 




Figure 2: Triangle and bubble topologies for the six-point. 


The finite colour-ordered amplitude takes the form 

(1^,2^, 3^, 4^, 5^, 6^) = Ci 23|4|5|6;4 fl23|4|5|6 
+ Cl2|34|5|6;4 fl2|34|5|6 ^^12|3|45|6;4 fl2|3|45|6 

-H cyclic perms. (33) 

Which agrees with 



Cl|2|3|H;0 - --2l4‘^|2(l^’2^,3^,H)si2S23 > 

Cl|2|3|H;4 = 0, (37a) 


As done for the five-point, we also consider other 
helicity configurations, A6(l“,2+,3^,4+,5^,6^), 
Afi (1-, 2-, 3+, 4+, 5+, 6+), As (r, 2-, 3', 4+, 5+, 6+), 
As (1-, 2+, 3 -, 4+, 5 -, 6+), As (r, 2+, 3', 4+, 5+, 6+), 
As (r, 2+, 3+, 4 -, 5+, 6+) and As (r, 2+, 3+,4”, 5', 6+), 
where we also have to consider contributions from 
triangles and bubbles, such topologies are depicted in 
fig|2] Moreover, the relevant tree-level are the ones that 
appear in eq. (l30l l and the six-points. 


GG^gggg, SS^gggg, (34) 


where the numerical value of each coefficient in eq. (l20l i 
agrees with NJet for those helicity configurations and 
with previous results[19, ^ 26, 2^371. 

5.3. The gggH Amplitude 


As final example we show the calculation of the 
leading colour-ordered one-loop helicity amplitude 
As (1^, 2^, 3^, H) in the heavy top mass limit. Since this 
amplitude is symmetric under cyclic and non-cyclic per¬ 
mutations of the particles we only consider the indepen¬ 
dent topologies for boxes, triangles and bubbles. 

The leading-order contribution reads as follows 


^4,H 


(l'",2+,3+,H) = 


<12><23><31)’ 


(36) 


The triple cut with two massive channels is 



H H 


Cl2|3|H;0 = -214 _h(1^>2^,3^,H)(5i3 + 533 ) , 

Cl2|3|//;2 = 0 , (38a) 

while the one with one massive channel only reads as 
follows: 



1 “ 1 “ 


Cl|2|3H;0 - 0 , 

CI|2|3H;0 = -2Af^ (1^ 2+, 3^ H) , (39a) 


Finally, the double cut is given by: 


C\2\iH - 




The quadruple cut is given by: 
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Cl2|3H;0 - 0 , 

C12|3H;2 (40a) 

The cut Ci 23 |H does not give any contribution. The re¬ 
quired trees to compute these coefficients are (l30l l and 


There are many outlooks for this job; they involve the 
computation of the analytical expressions for Higgs -H 
2(3) jets in the final state, and also the two-loop imple¬ 
mentation of our formalism. 


Acknowledgements 


GG^H, SS^H, 

GG^gH, SS^gH (41) 


where the Feynman rules for the Higgs-gluon and 
Higgs-scalar couplings in the FDF are given in Ap¬ 
pendix C of 101 . 

Then the colour-ordered one-loop helicity amplitude 
takes the form 


^4H°'’(1^>2^,3^,H) - Ci|2|3|H;0fl|2|3|H;0 


+ Ci2|3|H;0fl2|3|H;0 + Ci|2|3H;2 fl|2|3H;2 

+ ci 2 | 3 H ;2 /i 2 | 3 H ;2 + cyclic perms. (42) 


Which agrees with OSlI . 


The procedure for computing the one-loop ampli¬ 
tudes given above has been fully automated. In partic¬ 
ular, we have implemented the FDF Feynman rules (in¬ 
cluding the (-2e)-SRs) in FeynArts/FeynCalc 1391], to 
build automatically the tree-level amplitudes to be sewn 
in the cuts. Then, the coefficients of the master integrals 
are determined by applying the integrand reduction via 
Laurent expansion l40l] . which has been im plem ented in 
Mathematica, by using the package S@M 14111 . 


6. Conclusions 

At one-loop level, we have explored the unitarity 
methods and we have provided a new formalism with 
extended helicity spinors and consequently extended 
polarisation vectors, which allows for fully reconstruct¬ 
ing the full one loop scattering amplitude. We remark 
that there is an unified formalism in which the cut- 
constructible part and the rational part of a scattering 
amplitude can be found at once. It is enough just to 
give off-shellness to the internal momentum in a natural 
way related to the dimensional regularisation and then 
to perform multiple unitarity cuts for massive internal 
legs, where for a massless theory such a mass is exactly 
the off-shellness. 

We have presented a set of very non-trivial examples, 
showing that FDF scheme is suitable for computing im¬ 
portant 2 —> 2,3,4 partonic amplitudes at the next-to- 
leading order. 
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